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Abstract. We review the recent progress in studying the anomalous electroweak quartic 
gauge boson couplings (QGBCs) at the LHC and the next generation high energy e ± e~ 
linear colliders (LCs). The main focus is put onto the strong electroweak symmetry 
breaking scenario in which the non-decoupling guarantees sizable new physics effects 
for the QGBCs. After commenting upon the current low energy indirect bounds and 
summarizing the theoretical patterns of QGBCs predicted by the typical resonance/non- 
resonance models, we review our systematic model-independent analysis on bounding 
them via W^VF-fusion and WW Z / Z Z Z-production. The interplay of the two production 
mechanisms and the important role of the beam-polarization at the LCs are emphasized. 
The same physics may be similarly and better studied at a multi-TeV muon collider with 
high luminosity. 

PACS number(s): 11.30.Qc, 11.15.Ex, 12.15.Ji, 14.70.Pw 



1. Introduction 

The non-Abelian gauge structure of the standard model (SM) predicts the pres- 
ence of electroweak quartic gauge boson couplings (QGBCs) besides the couplings of 
triple gauge bosons. The electroweak symmetry breaking (EWSB) sector involves the 
would-be Goldstone boson [1] dynamics which generates the longitudinal components 
for W^^sothat they acquire the observed masses. Despite the astonishing success 
of the SM at scales up to O(lOOGeV) [2,3], this EWSB sector remains unverified [4]. 
Any new physics in the underlying Goldstone boson dynamics will cause the gauge 
boson self-interactions to deviate from the SM. The quartic gauge boson interactions 
are particularly interesting because they can involve four longitudinal components 
which, according to the equivalence theorem [5,6], manifest at high energies as pure 
Goldstone boson interactions (that is independent of the SM gauge couplings). To 
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End of a Muon Collider " , November 6-9, 1997, at Fermi National Accelerator Laboratory, Batavia, 
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unambiguously test the couplings of the quartic gauge boson interactions (QGBCs), 
the high energy H^H^-fusion and triple gauge boson production processes have to be 
used, where the QGBCs directly appear at the tree level. It is therefore important 
to study how the future high energy colliders (such as the CERN LHC and e ± e~ 
linear colliders [7,8]) can sensitively probe the QGBCs for unveiling the mystery of 
the EWSB mechanism. 

The EWSB sector can interact weakly or strongly. The weakly coupled case (such 
as supersymmetric models [9]) ensures the new physics at higher scales to have de- 
coupling property [10] at low scales, while in the strongly interacting scenario [11] the 
nondecoupling guarantees the new physics scale to lie below or at 4nv ~ 3 TeV [12]. 
In the former case the light Higgs boson(s) plus superpartners have to be first dis- 
covered, while for the latter we expect sizable new physics deviations showing up 
in the quartic (and triple) gauge boson couplings, which is the focus of this review. 
Below the new physics scale A, all the new physics effects in the EWSB sector can be 
parametrized by a complete set of the next-to-leading order (NLO) effective opera- 
tors of the electroweak chiral Lagrangian (EWCL) [13], in which the SU(2)l®U(1)y 
gauge symmetry is nonlinearly realized 1 . Without experimental observation on any 
new light resonance [2,3], this effective field theory approach [14,12] provides the 
most economic description of the possible new physics effects. Among the complete 
set of the fifteen NLO operators, five of them characterize only the quartic gauge 
interactions [13]: 

[ U = Mi) 2 [Tr(V M V,)] 2 , \ 

. A = 4(i) 2 [Tr(V M V^)] 2 ; 

' C 6 = 4 (f) Vr(V M K)]Tr(TV^)Tr(TV^) , ' 

< C 7 = £ 7 (f) 2 [Tr(V^)]Tr(TV,)Tr(TV^) , > 

k Ao = ^io (x) 2 |[Tr(TV M )Tr(TV iy )] 2 . 

In (1), V„ = (D^U)W , D^U = d^U + W^U-UB^ , W M = igW^/2, B M = ig'B^/2, 
U = exp[iT a ir a /v] (with 7r a the would-be Goldstone boson field), and T = Ut^U^ is 
the custodial SU(2) C - violation operator. Here, the operators £ 4i5 conserve SU(2) C 
while £6,7,10 violate SU(2) C . The dependence on v and A is factorized out so 
that the dimensionless coefficient £ n of the operator C n is naturally of 0(1) [12]. 
Because they contain only QGBCs these five operators cannot be directly tested via 
their tree-level contributions at low energies and are therefore least constrained from 
the current data. So far, only some rough estimates have been made by inserting 
them into the one-loop corrections and keeping the log-terms only 2 . Here is a recent 

^ It is advised that whenever the decoupling theorem [10] becomes ineffective, the nonlinear real- 
ization should better apply. 

2 ) The ignored constant contributions plus the new loop counter-terms are of the same order of 
magnitude as the log-terms. So, some uncertainties (like a factor of 2 to 3) may naturally exist in 
these estimates. 
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estimate at 90% C.L. by choosing A = 2 TeV and setting only one parameter nonzero 
at a time [30,15]: 



-4 < U < 20 , - 10 < 4 < 50 ; , , 

-0.7 < 4 < 4 , - 5 < £ 7 < 26 , - 0.7 < £ 10 < 3 . ^ ' 

(2) shows that the bounds on the SU(2) C symmetric parameters £4,5 are about an 
order of magnitude above their natural size of 0(1); while the allowed range for the 
S'?7(2) c -breaking parameters 4-10 is about a factor of 0(10 — 100) larger than that for 
£ = ^Ap (= |JaT) derived from the p (or T) parameter: 0.052 < £ < 0.12 [15], 
for the same A and confidence level. To directly test the EWSB dynamics, it is 
crucial to probe these QGBCs at future high energy scattering processes where their 
contributions can be greatly enhanced due to the sensitive power-dependence on the 
scattering energy [15]. 



2. Quartic Gauge Boson Interactions and Underlying Models 



So far the full theory underlying this effective EWCL is not determined, it is thus 
important to analyze how the typical resonance/non- resonance models contribute to 
these EWSB parameters. Knowing the theoretical sizes and patterns of these param- 
eters tells how to use the phenomenological bounds derived in following sections for 
discriminating different new physics models. We mainly focus on the quartic gauge 
boson interactions (1) and consider [17] typical models such as a heavy scalar (S), a 
vector (V®) and an axial vector {A a ^} for the resonance scenario, and the new heavy 
doublet fermions for the non-resonance scenario. 



• A Non-SM Singlet Scalar Up to dimension-4 and including both SU(2) C 
conserving and breaking effects, we write down the most general Lagrangian for a 
singlet scalar which is invariant under the SM gauge group SU(2) L <g> U(1) Y - 



^eff 



d^Sd^S - M 2 S S 2 



V(S) 



2 4 



Tr [V^] 



'-^vS + '-^S 2 
2 4 



[TrTV M ] 



(3) 



where V(S) only contains Higgs self- interactions. The SM Higgs boson corresponds 
to a special parameter choice: k s — k' s — 1, k s = k' s = and V(S) = V(S)sm ■ A 
heavy scalar can be integrated out from low energy spectrum and the corresponding 
contributions to (1) are derived as: 



£1 = 0, £ s 5 = ^>0; €S = 0, £ s 7 = ^ 



(4) 



In (4), the deviation from k s = 1 and k s = signals a non-SM Higgs boson. 
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• Vector and Axial- Vector Bosons The S'-parameter measurement at LEP 
disfavors the naive QCD-like dynamics for the EWSB [11], where the vector p TC is 
the lowest new resonance in the TeV regime. This suggests a necessity of including 
the axial-vector boson [18] in a general formalism for modeling the non-QCD-like 
dynamics. We consider the vector V® and axial-vector A a ^ fields as the weak isospin 
triplets of custodial SU (2) c . {V, A} transform under the SM global SU (2) c as =>■ 
% = KV^t A^^ = KA^t, where % = V*r a /2, \ = A*T a /2, and E„ G 
SU(2) C . If {V, A} are further regarded as gauge fields of a new local hidden symmetry 
group H = SU(2)' L <g) SU(2)' R (with a discrete left-right parity) [18], we can write down 
the following general Lagrangian (up to two derivatives), in the unitary gauge of the 
group 7i 3 and with both S77(2) c -conserving and -breaking effects included 4 , 



rVA _ rVA 
^"eff ~~ "^kinetic 



K Q TxVl + KiTr (J* - 2V ll f + k 2 Tv (j£ + 2A„) 2 + k^A* 



TrTV, 



+ K\ 



2V a 



+ K 2 



+ K 3 



TiTA 



(5) 



where 



J v 



J L 



J,. = J, z 



and, by definition, = igV^ 
U^V^U, f = r 3 = C/tTC7, and 
7i. Among the above four new 
by normalizing the Goldstone 
eliminating the V and A fields 



= ^ f = e + 5^ f ) 

= igV«r a /2, Ap = igA^ = i~gA«T a /2, V„ = U^D^U = 
U = £ 2 . Here g is the gauge coupling of the group 
S'f/(2) c -conserving parameters « n 's, /to is determined 
kinematic term: k = — 4ft 2 «; 3 /(4fi; 2 + ^3) • After 
in the heavy mass expansion, we derive £ n 's below: 





' h = 


^1 + ^2 


' w — 
£4 — 


/W 


= l/pv^sfvA- 1 ] 2 > 




h = 






fa 


= [r/ 2 (r/ 2 - 2) + 16r/ 2 ] /fiV^gvA- 1 } 2 


< 


k = 




• — 


pv _ 
7 — 







£7 = 


e v 7 + £% 


/7a _ 
£ 6 — 




= -fj [4(3 - r/ 2 )r] + (1 - g 2 )v] /[2\/2gvK 






_ pv 1 pa 


pv _ 


- no 


= 


in 


which 











(6) 



7/ 



4k 2 



4k 2 + ac 3 



7/ 



2k 2 + AR 2 



2k 2 



(Ak 2 + k 3 ) + 2 (4k 2 + « 3 ) 4k 2 + /t 3 



(7) 



3 ) By "unitary gauge" we mean a gauge containing no new Goldstone boson other than the three 
ones for generating the longitudinal components of the known W, Z. In fact, it is not essentially 
necessary to introduce such a new local symmetry H for {V, A} [19] since H has to be broken 
anyway and {V,A} can be treated as matter fields [20]. The hidden local symmetry formalism is 
more restrictive on the allowed free-parameters (k„'s etc) due to the additional assumption about 
that new local group H. 

4 ) In the literature [18], only the S[/(2) c -conserving operators were given. 
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and A = min{My, M^} ■ After ignoring the SM gauge couplings g and g', 
{My, M^} ~ [gv^/Ki, gv^K2 + K3/4}, at the leading order. In (6), the factor 
1/lgvA^ 1 } 2 ~ Ki(A/My) 2 = O(ki) and all >S?7(2) c -breaking terms depend on fj . 
Note that the S'?7(2) c -symmetric contribution from the axial-vector boson interac- 
tions to 1% = —£% becomes negative for \r\\ < y/2 , while the summed contribution 
£ 4 = -£ 5 = [{jf - l) 2 + 16r) 2 ] /fiy/lgvA- 1 ] 2 > . The deviation of r] and/or fj from 
77(77) = represents the non-QCD-like EWSB dynamics. 



• Heavy Doublet Fermions Take for instance a model of one flavor heavy 
chiral fermions which form a left-handed weak doublet (Ul,D l ) t and right-handed 
singlets {Ur,D r }, and joins a new strong SU(N) gauge group in its fundamental 
representation. Their small mass-splitting breaks the SU(2) C and is characterized by 
the parameter cu — 1 — [Mu/Mof ■ The anomaly-cancellation is ensured by assign- 
ing the {U,D} electric charges as {+§, — \} ■ By taking {U,D} as the source of the 
EWSB, the W, Z masses can be generated by heavy fermion loops. The new contri- 
butions to the quartic gauge couplings of W/Z come from the non-resonant {U,D} 
box-diagrams. The leading results in the l/M UjD and uj expansions are summarized 
below: 



-2£l 



A 



1 2 



4:71V 



N 
12 



>0 ; 



A 



Attv 



7N 2 
240^ ' 



4o = ; (8) 



in which A = minjMjy, M D } . 



3. A Global Analysis on Probing QGBCs versus TGBCs at the LHC 



The general EWCL formalism [13] contains in total 15 NLO new operators whose 
coefficients (^ n 's) depend on the details of the underlying dynamics as exemplified in 
the previous section. It is shown [15] that, except for £0,1,8 (S,T, U), the current data 
only bound a few triple gauge boson couplings (TGBCs) to O(10) at the lcr-level 
and give no direct tree-level bound on QGBCs. The rough estimates of the bounds 
from 1-loop corrections still allow QGBCs to be of 0(5 — 50). For a complete test 
of the EWSB sector in discriminating different dynamical models, all these TGBCs 
and QGBCs (£„'s) have to be measured through various high energy VV-fusion and 
//^-annihilation processes. (V a = W ± , Z°.) For this purpose, a systematic global 
analysis [15] has been carried out which reveals the important overall physical pictures 
and guides us for further elaborate precise numerical studies (cf. Sees. 4-5). In 
performing such a global analysis we developed a precise electroweak power counting 
rule (a la Weinberg) for conveniently estimating all high energy scattering amplitudes 
and formulated the equivalence theorem (ET) [6] as a necessary physical criterion for 
sensitively probing the EWSB dynamics. Applying this counting method, we have 
carried out a systematic analysis for all V a V b -> V c V d and //W -> V a V h , V a V b V c 
processes by estimating the contributions to their ^-matrix elements from both the 
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leading order operator up to one-loop and the other 15 NLO operators at the tree- 
level. Based upon the basic features of the chiral perturbation expansion, we further 
build the following electroweak power counting hierarchy for the S'-matrix elements, 5 

3 E U 
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In the typical TeV region, for E G (750 GeV, 1.5 TeV), this gives: 

(9.3, 37) > [(0.55, 8.8), (2.0, 4.0)] > [(0.12, 0.93), (0.42, 0.42)] > 
[(0.025, 0.099), (0.089, 0.045)] > [(5.3, 10.5), (19.0, 4.7)] x 10^ 3 > (1.1, 1.1) x 10~ 3 , 

where E is taken to be the invariant mass of the VV pair. This power counting 
hierarchy can be nicely understood. In (9), from left to right, the hierarchy is built 
up by increasing either the number of derivatives (i.e. power of E/A) or the number 
of external transverse gauge boson Vr's (i.e. the power of gauge couplings). This 
power counting hierarchy provides us a theoretical base to classify all the relevant 
scattering amplitudes in terms of the three essential parameters E, f n and A plus 
possible gauge/ Yukawa coupling constants. 

At the event-rate-level, we have adopted the usual leading-log effective vector boson 
method [24] to reasonably and conveniently estimate the ^^-luminosities. In Fig. 1, 
the rate \Rb\ denotes an intrinsic background defined via the formulation of the ET 
as a necessary criterion for the sensitivities to the EWSB [6,15]. Fig. 1 shows that, 
at the 14TeV LHC with JC = lOOfb- 1 Luminosity and for A = 2TeV, the W + W + - 
fusion is most sensitive to £4,5 (QGBCs) and marginally sensitive to £3,9,11,12; while the 
qq' — > W + Z annihilation can best probe £3,11,12 and marginally test ^8,9,14- Hence, the 
W-fusions and //^-annihilations are complementary in probing the different sets of 
these NLO parameters (the QGBCs and TGBCs) at the LHC. 



4. Measuring the QGBCs via VFH^-Fusion Processes 

Though the LHC will give the first direct test on these new quartic gauge boson 
couplings (QGBCs), the large backgrounds limit its sensitivity and cutting off the 
backgrounds significantly reduces the event rate. As shown in Ref. [21], for the 
non-resonance W ± W ± production channels in the TeV regime only around 10 signal 
events were predicted at the LHC with a 100 fb _1 annual luminosity after imposing 
necessary cuts in the gold-plated modes (by pure leptonic decays). The corresponding 
study at the TeV e ± e~ LCs opens a more exciting possibility [22,23]. 

In this and next sections we review how to make further precision constraints 
for the QGBCs via the WW-fusion [25,26], WWZ/ZZZ-pro duct ion [17] 6 , and their 

5 ) For //W — > VVV amplitudes, there is an additional factor 1//^ by dimentional counting. 

6 ) The WWZ/ZZZ-production in the SM was first studied in Ref. [27], and later some analyses 
on including the anomalous couplings have also appeared [28] for the case of unpolarized e T beams. 
For a very recent study similar to Ref. [17] for WWZ/ZZZ-pyoductioia, see Ref. [30]. 
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interplay at LCs [17,29], which is much cleaner than the LHC so that the final state 
W/Z's can be detected via the dijet mode and with large branching ratios. Due to the 
limited calorimeter energy resolution, the misidentification probability of W versus 
Z and the rejection of certain fraction of diboson events should be considered [23]. 
Inclusion of the leptonic decay of Z to e~e + and is also useful. According to 

the study of Ref. [23], the (mis) identification probabilities of W and Z via jet-decay 
mode can be derived as 

W 85% W, 10% Z, 5% reject , 

(10) 

Z 22% W, 74% Z, 4% reject . v ; 

The detection efficiencies for WW, ZZ and WZ final states are thus estimated below, 
which are about 34%: 

e ww = [0.68 x 0.85] 2 = 33.4% , e ZZ = [0.70 x 0.74 + 0.067] 2 = 34.2% , 

e wz = [0.68 x 0.85] [0.70 x 0.74 + 0.067] = 33.8% . ^ ' 

To completely determine all the QGBCs, we need at least five independent pro- 
cesses. From W^W^-fusions alone, we can have 

Full process : Sub — process : Relevant parameter : 

e-e+ -> vuW~W + , (W~W + -> W~W + ), (£4,5) , 

e - e - _> WW-W- , {W~W~ -> W~W~), (£4,5) ; ^ 

e - e + ^ V vZZ , (W~W+ -> ZZ), (£ 4 , 5 ; 4,r) , 

e - e + e ±^T Z ? ( W t Z ^ ^t Z)) (^ 5; ? 

e-e+^e-e+ZZ, (ZZ -> ZZ), ([£ A + £ 5 ] + 2[£ 6 + £ 7 + i w \) ; 

where in the round backets the corresponding fusion (signal) sub-processes are given. 
We see that £ 4j5 can be cleanly tested via the first two processes in (12), as shown 
by Fig. 2. (All our plots have chosen the new physics cutoff as A = 2 TeV and 
the numerical results for other values of A can be obtained via re-scaling.) By 
including the third and fourth reactions £qj can be further disentangled. Finally 
the fifth channel provides a unique probe on £ w . Though this scheme is complete 
in principle, the realistic situation is more involved. Note that the rate of the last 
reaction in (12) is significantly lower than all others due to the double suppressions of 
the e-e-Z couplings while the fourth channel has huge backgrounds which are uneasy 
to overcome [23,26]. But ZZ — > ZZ also has an advantage due to the absence fusion- 
type backgrounds and the triple gauge boson couplings have no contribution either. 
This makes it relatively cleaner than others. Since the parameter £10 appears only 
in 4Z vertex, the above last channel has to be used anyway when only the fusion 
mechanism is studied. (For the process e~e + — > ZZZ on £ w , see Sec. 5.) Since the 
large backgrounds make the W^Z-channel less useful (see Fig. 4a below), we propose 
to use the production e~e + — ► WW Z (cf. Sec. 5) to complete this five parameter 
determination. 
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From the above analysis, we finally summarize below the 90% CL. (one-parameter) 
fusion-bounds for A = 2 TeV at a later stage of the LC with the energy ^/s =1.6 TeV 
and the integrated luminosity J C = 200 fb _1 : 

-0.13 < l A < 0.10 , - 0.08 < l 5 < 0.06 ; , . 

-0.22 < 4 < 0.22 , - 0.12 < l 7 < 0.10 , - 0.21 < £ w < 0.21 ; ^ 

which are very stringent. Here we have used a 90% (65%) polarization for the e"(e + ) 
beam. 



5. WW Z/ZZZ-Production and its Interplay with IW-Fusion 

To probe the QGBCs (1), we know [15] that the H^H^-fusion amplitudes have the 
highest E-power dependence in the TeV regime while the s-channel signals of the 
WWZ/ZZZ-pro duct ion lose an enhancement factor of (E/v) 2 relative to that of the 
fusion processes. When the collider energy is reduced by half (from 1.6 TeV down 
to 800 GeV), the sensitivity of the H^H^-fusion decreases by about a factor of 20 or 
more [25,26]. We thus expect that ee — > WWZ, ZZZ become more important at 
the earlier phase of the LCs and will be competitive with and complementary to the 
fusion processes for the later stages of the LCs around 0.8 ~ 1 TeV [17]. In fact, 
it was revealed that even at the 1.5/1.6 TeV, e + e~ — > WWZ plays a crucial role in 
achieving a clean five-parameter analysis [17,29]. 

To avoid the potential fusion backgrounds from e"e + — > eeZZ, eeWW, we now only 
add the Z — > decay besides the dijet-decay mode. The detection efficiencies 

for ZZZ and WWZ final states are thus estimated to be about 16.8% and 18.4%, 
respectively. It turns out that e~e + — > WWZ has huge backgrounds due to the t- 
channel v e or e-v e exchange, and the kinematic cuts alone help very little. However, we 
find that such type of backgrounds involve the left-handed W-e-u coupling and thus 
can be effectively suppressed by using the right (left) -hand polarized e~(e + ) beam. 
The highest sensitivity is reached by maximally polarizing both e~ and e + beams. The 
crucial roles of the beam polarization and the higher collider energy for the WW Z- 
production are demonstrated in Fig. 3a, where ±l<r exclusion contours for £4-^5 are 
displayed at y/s = 0.5, 0.8, 1.0 and 1.6 TeV, respectively. The beam polarization 
has much less impact on the ZZZ mode, due to the almost axial-vector type e-Z-e 
coupling. Including the same polarizations as in the case of the WWZ mode, we find 
about 10 — 20% improvements on the bounds from the ZZZ-pro duct ion. Assuming 
the two beam polarizations (90% e~ and 65% e + ), we summarize the final ±lcr 
bounds for both ZZZ and WWZ channels and their combined 90% C.L. contours 
for 0.5 TeV with / C = 50 fb -1 in Fig. 3b (representing the first direct probe at the 
LC) and for 1.6 TeV with J C = 200 fb _1 in Fig. 3c (representing the best sensitivity 
gained from the final stage of the LC with energy around 1.5/1.6 TeV). Note that, 
the 90% C.L. level bounds on i 4 -£ 5 at 0.5 TeV are within 0(10-20), while at 1.6 TeV 
they sensitively reach 0(1). The WWZ channel gives the same bounds for £4-^5 and 
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£§-£7, while the ZZZ channel imposes stronger bound on £§-£7 due to a factor of 2 
enhancement from the 4Z-vertex. £ w only contributes to ZZZ final state and can be 
probed at the similar level. 

For comparison, a parallel analysis to Fig. 3b-c is further performed for the case 
without e + -beam polarization (but with e~ polarization the same as before). For a 
two-parameter (£4,5) study, the results are listed below at 90% C.L.: 

at0.5TeV: -12 (-18) < £4 < 21 (27), - 17 (-22) < £ 5 < 9.5 (15); 

at 1.6 TeV : -0.50 (-0.67) < £4 < 1.5 (1.7), - 1.3 (-1.5) < £ 5 < 0.36 (0.58); 

(14) 

where the numbers in the parentheses denote the bounds from polarizing the e~- 
beam alone. The comparison in (14) shows that without e + -beam polarization, the 
sensitivity will decrease by about 15% — 60%. Therefore, making use of the possible 
e + -beam polarization with a degree around 65% is clearly helpful. In the above, 
the total rates are used to derive the numerical bounds. We have further studied 
the possible improvements by including different characteristic distributions, but no 
significant increase of the sensitivity is found. 

Now, we are ready to analyze the interplay with VLW-fusion processes. As noted in 
Sec. 4, the VFZ-channel in (12) has large 7-induced eeWW background in which one 
e is lost in the beam-pipe and one W misidentified as Z. A cut on the missing p±{v) 
is imposed to specially suppress this background. Even though, the final sensitivity 
still turns out to be less useful in constraining the £ e -£ 7 plane (cf. Fig. 4a) [26]. 
To sensitively bound {^6^7} (especially £ 6 ) well below 0(1), we propose to use the 
production e~e + — > WW Z. Fig. 4a demonstrates the interplay of VLW-fusion and 
WWZ-pro duct ion for discriminating the S'[/(2) c -breaking QGBCs £§-£7 at s/s = 
1.6 TeV. The ZZZ-production can also bound £ w , in addition to the eeZZ fusion- 
channel in (12). Assuming that £4^-6,7 are constrained by the processes mentioned 
above, we set their values to the reference point (zero) for simplicity and define the 
statistic significance S = \J\f — Aol/V^o which is a function of £ w . Here Af is 
the total event-number while Ao is the number at £10 = 0. As shown in Fig. 4b, 
at 1.6 TeV, the sensitivity of e~e + — > eeZZ for probing £i is better than that of 
e~e+ -> ZZZ. 

In summary, the first direct probe on these QGBCs will come from the early phase 
of the LC at 500 GeV, where the PFW-fusion processes are not useful. The two mech- 
anisms become more competitive and complementary at energies ^/s ~ 0.8 — 1 TeV. 
From (13) and Table 1, we see that at a later stage of the LC with ^/s = 1.6 TeV, 
the 90% C.L. one-parameter bounds on £4^ from WWZ/ZZZ-modes are about a 
factor of 3 ~ 6 weaker than that from VFVF-fusions; while the bounds on ^6,7,10 are 
comparable. In a complete multi-parameter analysis, the PFiyZ-channel is crucial 
for determining £ 6 -£ 7 even at a 1.6 TeV LC (cf. Fig. 4a). 



Table 1: Combined 90% C.L. bounds on £ 4 _i from WWZ/ZZZ-pmduction. For sim- 
plicity, we set one parameter to be nonzero at a time. The bound on £10 comes from 
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ZZZ-channel alone. 



(TeV) 
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/^(fb- 1 ) 
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WW 7, 17,7,7, 

V V VV ZJ I ZJ ZJ ZJ 

Bounds 
(at 90%C.L.) 


- 9.5 <£ 4 < 11.7 

-5.0 < 4 < 5.8 
-5.0 < £ 7 < 5.7 
-4.3 < £ w < 5.2 


- 2.7 < £ 4 < 3.2 

— 3 1 < ft < 9 3 

-1.5 < 4 < 1.6 
-1.5 < ^ 7 < 1.5 

-i.4 < e 10 < i.4 


- 1.7 < 4 < 2.0 
_i n < f c < i i 

-L . c/ _^ -CJ3 _ J- • t: 

-0.95 < 4 < 1.0 
-0.95 < £ 7 < 0.92 
-0.83 < £i < 0.88 


- 0.50 < £ 4 < 0.58 

-0.28 < 4 < 0.28 
-0.28 < ^7 < 0.26 
-0.26 < £i < 0.26 


Range of \£ n \ 


< 0(4 ~ 10) 


< 0(1 ~ 3) 


< 0(0.8 ~ 2) 


< 0(0.3 ~ 0.6) 



6. Concluding Remarks 

Despite the constantly increasing evidence in supporting the Standard Model (SM) 
over the past 30 years, we particle physicists have been struggling in search for New 
Physics Beyond the SM so far [2,3]. Among the numerous ways for going beyond the 
SM, the Higgs boson hypothesis [31] stands out. The updated direct Higgs search 
at LEP [2] puts a 95%C.L. lower bound m H > 89.3 GeV. Due to the discrepancy 
between the precision Z-decay asymmetry measurement and the direct Higgs search 
limit, the combined 95%C.L. upper Higgs mass bound from the global fit has been 
shown to significantly increase toward the TeV regime [4]. However, the unitarity [32] 
and triviality [33] theoretically forbid the SM Higgs mass to go beyond the TeV scale, 
at which we are facing an exciting strong electroweak symmetry breaking (EWSB) 
dynamics. Below the new heavy resonance, we have to first probe the EWSB pa- 
rameters formulated by means of the electroweak chiral Lagrangian (EWCL), among 
which the quartic gauge boson interactions penetrate the pure Goldstone dynamics. 
After commenting upon the low energy indirect bounds and analyzing the different 
patterns of these quartic couplings predicted by the typical resonance/non-resonance 
models, we estimate the sensitivity of the LHC to probing these couplings, and then 
analyze the constraints on them via WW Z / Z Z Z-production and HW-fusion at the 
next generation e ± e~ linear colliders (LCs). The interplay of the two production 
mechanisms and the important role of the beam-polarization at the LCs are revealed 
and stressed. 

Finally, we remark that the same physics may be similarly and better studied 
at a multi-TeV muon collider (MTMC) (y/s ~ 3 - 4 TeV) with high luminosity 
(~ 500-1000 fb^Vyear) [34,35]. Due to the higher center mass energy of the MTMC, 
certain unitarization on the EWCL is needed for studying the H^iy-fusions. The 
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other muon collider options, like and are likely to be as easily achieved 

as [i~fi + mode. Furthermore, the large muon mass relative to the electron mass makes 
the initial state photon-radiation of the muon collider much less severe than that of 
the electron collider. The two drawbacks of a muon collider in comparison with an 
electron linear collider are [35]: (i) substantial beam polarization (> 50%) can be 
achieved only with a significant sacrifice in luminosity; (ii) the 77 and fi'-f options are 
probably not very feasible. 
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FIGURE 1. A classification of the contributions from all 15 next-to-leading order operators at 
14 TcV LHC (with 100 ft" 1 annual luminosity) for A = 2 TcV. 
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FIGURE 2. Determining the 5 , C/(2) c -symmetric parameters ^4-^5 at 1.6 TeV e~e + /e~e~ LCs. 
Here the ±lcr exclusion contours are displayed, (a), unpolarized case; (b). the case with 90%(65%) 
polarized e~(e + ) beam. Contributions from three types of resonance models (scalar, vector and 
axial-vector) to (At, £5) are shown by the thick solid lines. The different points on these solid 
lines correspond to different values of their couplings to the weak gauge bosons. Note that for 
axial- vector- type, it is also possible to have £4 + £5 = with €4 > 0, i.e., similar to the vector-type 
case. This makes the discrimination more involved. Big-star: from a scalar; black-triangle-down: 
from a vector; black-lozenge: from mixed contributions of a heavy scalar and vector. (Here we 
typically set these heavy resonances around 2 TeV.) 
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FIGURE 3. Probing £4-^5 via WW Z and ZZZ production processes. The roles of the polarization 
and the higher collider energy for e~e + — > WWZ are shown by the ± la exclusion contours in (a). 
The integrated luminosities used here are 50 ftr 1 (at 500 GeV), 100 ftr 1 (at 800 GeV) and 200 ft)" 1 
(at 1.0 and 1.6 TeV). In (b) and (c), the ± la contours are displayed for ZZZ /WW Z final states 
at y/s =0.5 and 1.6 TeV respectively, with two beam polarizations (90% e~ and 65% e + ); the thick 
solid lines present the combined bounds at 90% C.L. 



16 





FIGURE 4. Interplay of the WW-fusion and WWZ/ZZZ-pioduction for discriminating the 
5£/(2) c -breaking parameters £$-£7 and £10 at -,/s =1.6 TeV with J C —200 fb — 1 : (a). ± la ex- 
clusion contours for e~e + — > vvZZ, e + vW ~ Z / e~ vW + Z , and e~e + — > WW Z with polarizations 
(90% e~ and 65% e + ). (b). Statistic significance versus £\q for e _ e + — > ZZZ, e~e + ZZ (with 
unpolarized e T beams). 
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